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December 24, 2003 

Abstract. We show that the atmospheric and magnetic height variations are coupled in general MHS equilibria with gravity 
when isolated thin non-force-free flux tubes are present. In gas-dominated environments, as in stellar photospheres, flux tubes 
must expand rapidly with height to maintain pressure balance with the cool surroundings. But in magnetically dominated envi- 
ronments, as in stellar coronae, the large-scale background magnetic field determines the average spreading of embedded flux 
tubes, and rigidly held flux tubes require a specific surrounding atmosphere with a unique temperature profile for equilibrium. 
The solar static equilibrium atmosphere exhibits correct transition-region properties and the accepted base coronal temperature 
for the sun's main magnetic spherical harmonic. Steady flows contribute to the overall pressure, so equilibria with accelerated 
wind outflows are possible as well. Flux tubes reflect a mathematical degeneracy in the form of non-force-free fields, which 
leads to coupling in general equilibrium conditions. The equilibrium state characterizes the system average in usual circum- 
stances and dynamics tend to maintain the MHS atmosphere. Outflows are produced everywhere external to rigidly held flux 
tubes that refill a depleted or cool atmosphere to the equilibrium gas profile, heating the gas compressively. 
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1. Introduction 



In our observation of solar 'threads' ( iNovember & Koutchmvl 
I 1 

11996 . hereafter NK96) we puzzled over whether coronal flux 
tubes should expand with height to maintain pressure balance 
with the ambient gas-pressure decrease or follow the lines of 
the external magnetic field. After all, if the corona is indepen- 
dently heated, the gas and magnetic height variations should 
be generally different. The many observations of EUV and X- 
ray loops as well as the purely density-sensitive observations of 

white-light > 15 Mm 'voids' (.MacOueen et al 1974X > 2 Mm 

threa ds, and scintillation > 1 kilometer filamen tary microstruc- 
ture jColes fc HarmoJl978tlWoo et alJll995l) all suggest that 
the large-scale solar coronal magnetic field is interspersed with 
flux tubes with a great range of sizes. Figure shows dark 
and bright threads in a square region about a solar radius on a 
side over the west limb taken with the Canada-France-Hawaii 
Telescope (CFHT) at the unique total-eclipse opportunity on 
Mauna Kea on July 11, 1991. The dark and bright threads, 
which were essentially unchanged over the 4-minute eclipse 
duration, appear to be aligned and organized in arched and ra- 
dial surfaces, probably current sheets, overlapping in their pro- 
jection onto the plane of the sky in the line-of-sight view. 

We know that isolated thin flux tubes must expand with 
height in balance with the ambient gas-pressure decrease. The 
flux-tube gas pressure difference, external minus internal, is re- 
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Fig. 1. Unsharp-masked time-averaged white-light eclipse im- 
age processed by subtracting a 6 arcsec — 4300 km Gaussian 
convolution shown with contrast range A In / = 0.8%. Dark and 
bright threads are probably flux tubes in current sheets. 

lated to the magnetic pressure difference by the static equilib- 
rium relation 
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The Wilson depression in sunspots shows that at least large- 
scale photospheric flux tubes are nearly evacuated. Evacuated 
flux tubes exhibit the maximum magnetic flux concentration 
and field strength possible. In evacuated flux tubes = 
and the photospheric external field - 0, so the inter- 
nal field strength must follow the ambient gas-pressure de- 
crease Bii/f- - ^jipeif), giving flux tubes that expand with 
height to conserve their total internal magnetic flux. Since the 
scale height for peif) is small in the relatively cool photo- 
sphere and chromosphere, a dramatic flux-tube expansion oc- 
curs, which is the accepted exp lanation for the solar 'mag- 
netic canopy' dOiovan elli' 1982). Both voids and a range of 
sizes of dark threads were also found to have contrasts con- 
sistent wit h what is expected fo r fully evacuated cylindrical 
flux tubes (MacOueen et al. 1983; NK96). In the corona, where 
the external magnetic field is not zero, the magnetic flux in- 
side evacuated flux tubes must decrease with height to fol- 
low the total external ambient gas plus magnetic pressures, or 
Bi(r)^ = %7ipe{r) + Beif)'. This gives a small flux-tube expan- 
sion with height in the large gas-pressure scale height of the hot 
corona. A slight flux-tube expansion with height seems evident 
in some of the threads in Fig. [2 

We can argue too that the radial expansion of isolat ed thin 
flux tu bes must conform to the external magnetic field. IParkeJ 
ill972h shows that 

(B,-V)B,=0, (2) 

in MHS equilibria without gravity, where B^ represents a uni- 
form background magnetic field, and the field in flux tubes B, 
contains all of the spatial variations. The theorem is locally ap- 
plicable everywhere in perturbed MHS equilibria with gravity, 
i.e. Eq. (|2} with Bg allowed to vary slowly spatially, as we dis- 
cuss in Sect. 12 Thus thin flux tubes in a gravitational field must 
spread to follow the lines of the relatively uniform background 
magnetic field; the radial functional variations must be directly 
related, that is \Bi{r)\ oc \Be{r)\ along flux tubes. 

Therefore both ways of looking at the problem are essen- 
tially correct. The coronal flux-tube field must comply with the 
external field B,(r)^ oc Be{r)^, and the radial variation of the 
external gas pressure outside isolated thin evacuated coronal 
flux tubes is proportional to the magnetic pressure difference 
Peir) oc Bi{r)^ - Bi,{rf-. Hence, the external gas and mag- 
netic pressures must be coupled with Pe{r) oc B^r)^ outside 
and along isolated thin evacuated flux tubes. The background 
magnetic field is dominant in lower stellar coronae and must 
determine the expansion of embedded flux tubes. The expan- 
sion of evacuated flux tubes then defines the radial variation 
of the ambient gas pressure for equilibrium, the gas pressure 
being horizontally stratified in the external force-free or po- 
tential field. More generally, flux tubes of arbitrary base pres- 
sure but in thermal balance with the local ambient atmosphere 
satisfy pi{r) oc pe{r), giving the same external atmosphere 
Peir) oc Beif)^. Thus an atmosphere containing isolated thin 
non-force-free flux tubes in different forms of local balance 
must exhibit a temperature profile defined by the dominant 
background magnetic field for a perfect gas in hydrostatic equi- 
librium. 



The properties of flux tubes reflect a natural degeneracy 
in the form of non-force-free magnetic fields, and imply the 
existence of an underlying mathematical relationship between 
an equilibrium atmosphere and large-scale magnetic field. In 
Sect. 121 we show that the height variations of the gas pressure 
and magnetic field amplitude are directly coupled in the most- 
significant-order terms of the gravitational perturbation in gen- 
eral non-force-free MHS equilibria, with smaller field diver- 
gence and curl effects. In low-/? systems, with strong magnetic 
fields, the external gas pressure corresponds to a small differ- 
ence Pe{r) = (Bi(r)^ - Be{r)^)/{Sn), which allows large pres- 
sure deviations for a given relative flexibility of the large-scale 
magnetic field around individual flux tubes. Magnetic fields 
can exhibit and persist in significant deviations from equilib- 
rium, however the coupling Peir) oc B^ir)^ describes the gen- 
eral equilibrium condition for arbitrary /3. 

Historical MHS equilibrium solutions with gravity have 
also suggested that a general relation ship exis t s betw een the 
large magnetic and gas-pressure scales. Dungev ( 1953') derived 
plane-parallel MHS equilibria for a magnetic field with i ts de- 
generate axis directed perp endicular to gravity (see also: ILowI 
ll975tlCheng & Choe'l998V His solutions are represented in a 
modified Grad-Shafranov (GS) Equation that retains the gas- 
pressure radial scale-height variation as a separable factor in 
the pressure term. Dungey showed that when the magnetic 
field is also zero in its degenerate axis, the separable factor 
can be absorbed into the coordinates, giving the classical GS 
Equation in a remapped coordinate system. The new GS coor- 
dinates are defined by the gas-pressure scale height, meaning 
that magnetic fields in a prespecified atmosphere must be spa- 
tially distorted. Horizontal flux tubes exhibit substantial radial 
deformation in a solar corona of given co nstant temperature 
llZweihel & Hundhausenlfl98i lT,owlll992h . Depending upon 
their size compared to the gas-pressure scale height, flux tubes 
are either compressed or expanded in the gravitational direc- 
tion, and approximately undistorted only in a specific atmo- 
sphere. 

In Sect. |3l we examine equilibrium radial profiles of tem- 
perature and flux-tube diameter using an approximate solution 
to the equilibrium equations. The profiles reproduce the salient 
features of the solar atmosphere: the rapid expansion of flux 
tubes in the photosphere where the atmospheric gas pressure 
is dominant, a transition region of reasonable thickness located 
just above the height where the gas and magnetic energy den- 
sities are equal, and a corona of correct base temperature for 
the sun's main magnetic spherical harmonic. Steady flows do 
not alter the coupling, but may produce a Bernoulli pressure, 
which modifies the gas pressure, and leads to possible acceler- 
ated steady wind solutions. 

Although the dynamics of the system can be quite com- 
plicated, the equilibrium state characterizes the average i n 
many types of MHD systems, e.g. 'quasi-steady ' ilT,owlll98nl) . 
Hydrodynamics in such systems tend to maintain the MHS at- 
mosphere, as we discuss in Sect. |3 The possible importance 
of MHS atmospheres for other problems in physics and astro- 
physics beyond the formation of solar and stellar coronae and 
in more marginal non-quasi-steady physical conditions needs 
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to be considered, but discussion lies beyond the scope of this 
paper. 

2. MHS equilibria with gravity 

The equations for static equilibrium in a gravitational field are 
the MHS equation and Gauss' Law 



Vp + -Vr = —{VxB)xB, 

h An 

VB = 0. 



(3) 
(4) 



The static contribution to the pressure scale height h is directly 
related to the temperature for a perfect gas, h = plipgir)) - 
kT l{mg(r)). The variables have their usual meanings: B is the 
magnetic field vector in Gaussian units, p the gas pressure, p 
the density, T the temperature, m the mean particle mass, k 
Boltzman's constant, c the speed of light (used below); g{r) - 
GMJr^ is the gravitational acceleration, with G the universal 
gravitational constant and M, the stellar mass. The coordinate 
r measures the distance from the system or stellar center, and 
the formal notation Vr is adopted for algebraic convenience to 
denote the radially directed unit vector. 

The set of four vector-element Eqs. and (0} contain five 
unknowns, p, h, and the three vector elements of B. However, 
the solution space for B is more restrictive than would be ob- 
tained with linear relations, as we discuss in this Section and 
elaborate in AppendixlAl 

The gravitational term in Eq. Q compared to the pressure 
gradient is of order d/h, where is a characteristic flux-tube 
thickness; for the main power in threads in the solar corona 
dih < 10"^. In the vicinity of non-force-free fields gas-pressure 
changes must be mainly magnetically determined, and only far 
away can the gravitational gradient be significant. Where the 
Lorentz term is negligible, the gas pressure is hydrostatic and 
defined by a single radial temperature profile. Since the gravi- 
tational term is small in non-force-free conditions, the classical 
equilibrium solutions without gravity described in Appendices 
|A]and|B]are applicable, and Eqs. Q and (|2} are valid in ev- 
ery small volume, but with the allowance that quantities may 
vary slowly spatially. Thus we obtain coupled atmospheric and 
magnetic height variations, at least if certain types of flux tubes 
are present as described in the Introduction. 

Conditions along field lines are represented by a parallel- 
field equation 



B 



0, 



(5) 



which is the projected component of the MHS Eq. ^ along 
B. The relation says that gas pressure changes along field lines 
must be hydrostatic, which can be seen by introducing the in- 
tegrating factor (p into the gas pressure 

P = (pp, (6) 
for an unrestricted function p. Upon substitution 

B- (^(pVp + p iv</) + jVrJj = 0, (7) 



so taking for the integrating factor 



1. 



Vlni* = — Vr, 

h 



(8) 



we obtain B-Vp - 0, or p constant along field lines. The only 
solutions to Eq. ^ are the ID (f> - (p{r) and h - h(r) or T - 
T(r), represented in the hydrostatic relation 



din (pir) _ 1 
dr h{ry 

which can be integrated to give 



{r) - exp 



(9) 



(10) 



where r, denotes an arbitrary base height for the radial varia- 
tions. Thus pressure variations along arbitrarily directed field 
lines follow the ID hydrostatic scale-height function (p(r) from 
a constant base pressure p\ 3D effects are possible because the 
base pressure p and the scale-height function (p{r) or temper- 
ature profile T{r) can vary from field line to field line. One 
scale-height function ^{r) with temperature profile T{r) ap- 
plies throughout force-free regions and in isolated non-force- 
free flux tubes in thermal balance with their local surroundings. 

Substituting Eq. (|6} into the MHS Eq. (|3} and using Eq. (|8} 
gives the cross-field equation 



1 

(r)V/3 = — (VxB)xB, 
An 



(11) 



for p allowed to vary from field line to field line. We adopt 
the single scale-height function 0(r), assuming for our argu- 
ment that one temperature profile applies throughout. At a base 
height Vs, <p(rs) - 1, and the cross-field Eq. (^J resembles the 
classical MHS Eq. iA.l\ without gravity. Since the base height 
rs is arbitrary, it is again evident that the classical equilibrium 
solutions must be applicable in every small volume. 

Since the base pressure p is constant along field lines, the 
variation in the scale-height function 0(r) on the left side of 
Eq. il 1> must be reproduced on the right side in the Lorentz 
term, and so in B^ too. Thus B must vary like (f>{r)^^'^ along 
field lines in non-force-free fields within djh, which is the or- 
der of the derivatives of the scale-height function. Formally we 
can absorb (p(r) into the magnetic field on the right side by de- 
composing B into the general product 



B = cPBir)B = 4>B2iry'^B = 0(r)"^^(r)B 



(12) 



A ID magnetic scale-height function ^^(r) is sufficient to can- 
cel 0(r); we introduce the residual radial multiplier ^(r) to al- 
low differing scale-height functions, leaving B an unrestricted 
vector function. Substituting back into Eqs. il 1> and elimi- 
nates (p(r) and gives 



Vr 

2h(r) 



(13) 
(14) 



The equations can be solved by substituting perturbative 
expansions in powers of the small quantity d/h{r), p - 
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Y,n PnidlHr))" and B = E« B„(d/h(r))" for n > with d 
a constant characteristic flux-tube thickness. The expansions 
separate according to powers of {d/h(r))", giving the most- 
significant order-zero equations 

Vpo = ^ (Vx^(r)Bo) xar)Bo, (15) 
V.f(r)Bo=0. (16) 

Equations (I15> and (I16> in /3o and ^(r)Bo are equivalent 
to the classical plane-parallel MHS equations without gravity, 
which are discussed in AppendixlAl The solutions are the clas- 
sical equilibria, which exhibit a degenerate direction z in the 
local coordinate system (x, y, z) taken to be arbitrarily oriented 
with respect to the radial r. There is no loss in generality in just 
taking ^(r) = 1 and (pBiif) - <pir) in Eq. il2\ leaving Bo the un- 
perturbed solution from Eq. iA.5\ and constant in z. If z / r the 
solutions for B exhibit a variation along the field direction that 
goes like (p(ry^^. In the case considered by Dunaev ( 1953) with 
z ± r, distortion effects are introduced perpendicular to the 
main magnetic field direction as described in the Introduction. 
Thus liB{r) - 2h{r) at every r, where the magnetic scale height 
is defined in the usual way, hB{r) - -[51n0B(r)/5r]"^ The 
scale heights are coupled to most-significant order in the per- 
turbing gravitational field in all non-force-free MHS equilibria 
containing thin flux tubes in thermal balance with their local 
surroundings. 

Additional perturbative equations for the base functions p 
and B originate with the rescaling of the magnetic field vector 
in r as discussed in Appendix Icl The rescaling leads to the 
small term Vr/(2/i(r)) in both Eqs. ( I13> and MAI . The term 
appears on the right side of Eq. J14> and is needed since a slow 
radial decrease in |B| cannot occur without a divergence of field 
lines; it represents a compensating creation of flux in B{r) with 
r. The term is subtracted from the curl in the Lorentz force 
in Eq. M3\ and corresponds to a small residual pseudo-force 
directed between B and the outward radial Vr; it represents 
a progressive left-handed twist of field lines at the rate of 1 
radian in every magnetic scale height Ihif). In MHS equilibria 
with gravity, the amplitude, the divergence, and the curl of B 
all vary on the large scale 2h(r). 

In principle anyway, a single thin flux tube somewhere in 
the atmosphere, evacuated or in thermal balance with its sur- 
roundings, fixes the equilibrium gas pressure for the entire at- 
mosphere, but all non-force-free fields must be consistent to 
satisfy mutually the horizontal pressure boundary condition. In 
reality of course, magnetic fields may exhibit local deviations, 
and the equilibrium state might only be reflected in the average 
in an atmosphere containing many flux tubes. 

A flux tube in thermal balance satisfies the condition 
Pe{r) oc Bei/f whether of increased or decreased field strength 
compared to the local background magnetic field. A flux tube 
with increased field strength |B,(r)| > |Be(r)|, must exhibit a 
decreased internal gas pressure in Eq. {I), p,(r) < Peif) consis- 
tent with the extreme dark-thread evacuated case of pi{r) - 0, 
which requires the minimum surrounding atmosphere for a 
given magnetic pressure difference. Flux tubes with less field 
strength than the external magnetic field |B,(r)| < |Bf (r)|, must 



have an increased gas pressure /?,(r) > pdr), consistent with 
bright threads. 

Threads appear to be isolated thin ffux tubes aligned with 
the background magnetic field and embedded in current sheets 
in a bimodal amplitude distribution, consistent with the natural 
equilibrium form discussed in AppendixlE] The temperature of 
bright threads was indeed found to be the same as the external 
surroundings within observational uncertainties (NK96). The 
solar EUV coronal temperature for unresolved quie t and ac- 
tive regions is relatively uniform, 1.5-2.1 x 10'' K liWithbrod 
^975), and a narrow temperature range around 2 x 10^ K is 
found for a broad range of sizes and dens ities of X -ray loops for 
all but prominence flare loops jPavis et al...l975i: IVaiana et alJ 
Rosner et al. 1978). Large-area coronal temperatures de- 
rived from the solar white-light scale height and forbidden line 
ratios give temperatures somewhat lower than the EUV and X- 
ray temperatures, i n the range of 1 - 2 x 10^ K at the solar radial 
height r = 1.15ro dOuhadiakurta et alJl992i) . 

A slightly heated (or cooled) flux tube exhibits an increased 
(or decreased) internal pressure pi{r) with height and thereby 
must affect a decreased (or increased) internal field strength 
|B,(r)| or an increased (or decreased) external pressure Pe{r) to 
satisfy the equilibrium Eq. ([ij. Small changes in internal field 
strength might be accommodated by changes in the flux-tube 
size as a state of local stress in the large-scale magnetic field, 
and changes in the surrounding ambient pressure pdr) might 
be modified by compensating external dynamics like we de- 
scribe in Sect.|3] Hot flux tubes should tend to cool to the tem- 
perature of their immediate surroundings, bringing them back 
to the normal MHS equilibrium condition. A relatively cool 
nonevacuated flux tube with excess field strength tends to evac- 
uation with height as a different equilibrium condition, and so 
might just persist out of equilibrium within about one internal 
scale height of its base. 

Cont inuous 3D non-f orce-free equilibrium solutions are 
claimed ( Low 198?," 1991), but these allow an arbitrary hydro- 
static ID gas pressure independent of the magn etic field (s ee 
the use of Po(r) and discussion ar ound Eq. (24 ^ oflLowll99lL i n 
Sect. Ill in Bogdan & Low 1986, or Sect. 4 in"Neukirch'l997|). 
In the general force balance, that is the cross- field Eq. ^ with 
0(r) allowed to vary from field line to field line, the pressure 
and magnetic variations cannot be separated except in force- 
free regions unconstrained by a horizontal pressure boundary 
condition, and so necessarily lacking isolated thin non-force- 
free flux tubes evacuated or in thermal balance with their sur- 
roundings. 

We have confined our study here to the static solutions, but 
it is straightforward to broaden consideration to include certain 
types of steady flows. Steady uniform flows along field lines 
add a flow pressure to the parallel-field Eq. (|5}, which leads to 
a modified hydro-steady relation in Eqs. (|9j and MQ\ . but the 
same cross-field Eq. (II U results. Thus we obtain a gas pressure 
containing static and wind components coupled to the large- 
scale magnetic field outside flux tubes. For a given gas pressure 
variation or scale-height function 0(r) as is defined by the mag- 
netic field in a magnetically dominated environment, the static 
contribution to the pressure scale height hi/) and correspond- 
ing temperature kT{r) - mg{r)h{r) are always larger with a 



L. J. November: MHS atmospheres 



5 



Steady flow than without (| ParkeJl960h . Steady flux-tube flows 
can modify the coupling too. 

3. Model atmospheres 

The derivation of an explicit static equilibrium solution is a 
formidable task because such an analysis must consider the 
specific form that distorted magnetic fields take in the pres- 
ence of a gas. However an approximate coupled scale-height 
function can be written 



Ps 



(17) 



where 0fl(r) and (f>B2(f) are nominal separate gas and mag- 
netic pressure solutions to the MHS Eqs. Q and @ normal- 
ized at a base height rj; the coefficient l/yS, defines the relative 
magnetic-field strength. The approximate coupled scale-height 
function 0(r) is proportional to the total pressure p + B^/ (Stt) of 
the nominal separate solutions. The gradient of the total pres- 
sure approximates the gas-pressure gradient minus the Lorentz 
term in the MHS equilibrium Eq. reflecting the local bal- 
ance in Eq. ([0, as discussed at the end of AppendixlEI The cou- 
pled scale-height function ^(r) thus gives a correct MHS solu- 
tion where one pressure component is dominant, p » B^/(87r) 
or B^/(8n) » p. 

At a base height taken as the stellar photospheric surface, 
we assume that magnetic fields are relatively weak, /3s » 1 
and 0(r) ^ <pa(f)- There the gas pressure is dominant and must 
follow its nominal atmospheric form not being much affected 
by the magnetic field, whereas the magnetic field lines must 
be highly distorted. The nominal atmospheric-pressure scale- 
height function (pair) is defined by the hydrostatic relation Eq. 
JlOl l with a temperature found using a transfer equation that 
includes all of the usual energy input and loss mechanisms. 

The atmospheric contribution (pair) characteristically falls 
off much more rapidly than the magnetic pressure (pBiir), so 
around some transition height r, where (pain) - (pBiir^lPs, ^{r) 
changes from (paif) to <pB2{r)//3s over about one gas-pressure 
scale height. Above r, the magnetic field is dominant 0(r) ^ 
(pB2(r)lfis, and the field must follow its nominal form irrespec- 
tive of the gas pressure. There the gas pressure profile may be 
very distorted from its nominal form, representing an added 
atmospheric heating produced by MHS restoring flows as we 
discuss in Sect.|4] 

Both the magnetic and gas pressure variations must be dis- 
torted from their nominal separate forms around the intermedi- 
ate transition height r,. Preserving the total pressure p+B^/ (Stt) 
gives an approximate flux-conserving extrapolation below the 
magnetically dominated corona and an ostensible gas-pressure 
extrapolation above the photosphere and chromosphere. The 
total relative pressure around the intermediate transition height 
goes from Aiot{<pair) + 0B2('")/y6j) for a nominal superposition 
of independent gas and magnetic pressures (pair) and (pB2{r)//3s 
with Atot the total surface area, to (Atot - Avac)0('") + Ag^ir) for 
the coupled pressure 0(r) with Avac the evacuated area and Ag 
the magnetically filled area. The two pressure totals are equal in 
general only with 0(r) from Eq. ilH and when the field-filling 
regions are evacuated with Ag = Avac. Anyway the choice of 



coupled scale-height function is not too critical for our demon- 
stration, as smooth switching functions with correct asymptotic 
behavior in the photosphere and corona exhibit similar temper- 
ature profiles even around the transition height. 

For demonstration purposes, we take (pair) for a polytrope 
atmosphere as derived in Appendix ID! 



Mr) = 1 - 



r - 1 r, 

r hs 



(-7) 



(18) 



where = hir^) - kT J(mgs) denotes the surface scale height 
and F the ratio of specific heats or polytrope adiabat. 

A nominal scale-height function (pB2{r) is defined for a 
single-spherical-harmonic potential magnetic field. The poten- 
tial magnetic form seems to be consistent with what is observed 
in the lower solar corona fo r r < 1 .brQ ( Altschuler & Newkirk 
ll969HSchatten et al.lll969 ^. The potential magnetic field vec- 
tor is the gradient of a sum of scalar spherical-harmonic com- 
ponent functions, each of which is the separable product of a 
2D surface function and radial multiplier The resulting 

magnetic field in each spherical harmonic goes like 1/r^^^ in 
all its vector elements; a monopole £ - Q exhibits a ra- 
dial falloff, a dipole field ^ = 1, a falloff, etc. For a single 
spherical harmonic the magnetic energy density falls off like 



IJB2' 



(r. 



>=(7) ■ 



(19) 



With a superposition of spherical harmonics, a nonuniform 
radial dependence can occur, and a representative global- 
average spherical harmonic £ might be obtained by appropri- 
ately weighting flux-tube locations. 

The temperature is written from the hydrostatic Eq. (|5Jl 



kT(r) — mg(r)h(r) — —mg(r) 



1 d(p(r) 
<p(r) dr 



(20) 



Taking ^{r) to be the approximate ^(r) from Eq. Mil with (pair) 
from Eq. J18l l and (pB2{r) from Eq. il9\ . we obtain an explicit 
formula for T(r). In the corona (p{r) ^ (pB2{r)//3s, and the for- 
mula exhibits the limiting scale height 



and temperature 

GMsih mg(r)r 



kTcor(r) 



{2( + A)r 2( + A' 



(21) 



(22) 



The limiting base coronal temperature TcoArs) is defined en- 
tirely by the stellar surface gravity gs, radius r,, mean parti- 
cle mass m, and main magnetic spherical harmonic £. For a 
toroidal field £ - 2, the scale height at the stellar surface is 
hcm-irs) - Ts/S, which gives a solar base coronal temperature 
of T^-ojirs) - 1.73 X 10^ K (using ill = Q.6mp for nip the proton 
mass, m = 1 x 10^24 g, = = 6.96x10'° cm, g^ = 2.75 x lO" 
cm s"^, and k - 1.38 x 10"'^ ergs K~'). With a solar base 
coronal temperature of Tcor(rs) = 1.6 x 10^ K, we obtain the 
spherical harmonic ^ = 2.33. 

Figure |2 shows some example scale-height functions ^(r) 
from Eqs. ( fTTl - in the lowest part of the solar atmosphere 
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^^^^ ^ 1000 km ^ 

Fig. 2. Log of scale-height function logio ^C'') f"'" ^ solar poly- 
trope atmosphere with F = 1 . 1 and a potential magnetic field 
^ = 2.33 for different^,. 
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Fig. 3. Relative flux-tube diameter d{r)ld(rs) in lower solar at- 
mosphere with { =2.33 and =5000 for different polytrope 
adiabats F (upper panel), and with { =2.33 and F = 1.1 for 
different fis - 5, 50, . . . 5.E5 {lower panel). The arrow on each 
curve locates the most rapid temperature change for the model. 

around the transition height r,. Scale-height functions with the 
polytrope adiabat F = 1.1 and spherical harmonic { - 2.33 are 
plotted for different fis, using the solar parameters and surface 
photospheric temperature =682 K from the VAL model 
atmospheres dVernazza et aljfl98ll) . which defines the surface 
scale height - kTs/(mgs) for (pair) in Eq. ( I18> . As all of the 
curves are based on the same F, they coincide until the transi- 
tion height r, for the model is reached, and then switch rapidly 
to the nominal magnetic (pBiif), which appears to be relatively 
flat on the log scale. 

The radial expansion of a flux tube is another way to vi- 
sualize the field strength decrease in the atmosphere and the 
properties of the coupled scale-height function. For conserva- 
tion of the total flux through the cross-sectional area of a flux 
tube 7r(t//2)^|B|, the flux-tube diameter d{r) must increase with 
radial distance r, d(r) oc oc (f)(r)^^^'^. Figure |3l portrays 

the relative flux-tube diameter as a function of height d(r)/d(rs) 
for different solar atmospheric models denoted by F and /S,,- 

For a given polytrope adiabat F, all flux tubes exhibit ap- 
proximately the same relative shape up to a height that depends 
upon /3s, where the field lines straighten up. Stronger fields ex- 
ist with a lower transition height and exhibit less overall rel- 
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Fig. 4. Lower solar atmosphere logior(r) with t - 2.33 and 
fis =5000 for different polytrope adiabats F (upper panel), and 
with I - 2.33 and F = 1.1 for different ySj (lower panel), shown 
with the solar VAL A (open circles) and F (filled circles) atmo- 
spheres. 



ative expansion before straightening up. The canopy depends 
upon the surface distribution of the fields, but strong fields 
might give the appearance of a lower canopy height too consis- 
tent with observations of sunspots (Giovanelli & Jones 1983)- 
The straightening height is somewhat below the location of the 
fastest temperature change due to the differing dependencies: 
the field strength, which goes like 0(r)'^^, corrects from its at- 
mospherically determined form at r,, but the most rapid temper- 
ature change occurs significantly higher where the slope of the 
scale-height function flattens out, as T(r) oc -(51n0(r)/5r)"'. 

Figure 13 illustrates the temperature variation in the lower 
solar atmosphere using Eq. ( I20> with =6820 K from the 
VAL atmospheres. The upper panel shows the model atmo- 
spheres with t - 2.33 and fis =5000 for different polytrope 
adiabats F, and the lower panel the model atmospheres with 
f = 2.33 and F = 1 . 1 for different yS^. The VAL A and F model 
temperatures for relatively cool inner network and hot network 
bright points are also shown. 

In the solar photosphere, the nominal atmospheric poly- 
trope follows the VAL temperature roughly, giving the best 
tradeoff between photospheric and chromospheric tempera- 
tures with F ^ 1 . 1 . Of course the features of a real solar at- 
mosphere can never be well approximated by a polytrope. The 
shape of the temperature function through the chromosphere 
and transition height, where the scale-height function 0(r) goes 
from (f>a(r) to ^B'i(r)lfis, depends upon F and fis but not on 
I, which separately determines the base coronal temperature. 
Temperature profiles with fis - 5000 or 1000 give best agree- 
ment around the transition height with the VAL A or F atmo- 
spheres, respectively, averaging approximately through the Lyo' 
plateau. 

For flux-tube evacuation in the photosphere B(rs)^ - 
%jip(rs), and we obtain a maximum field strength of \B(rs)\ - 
1 7 1 5 G using p(rs) - 1.17x10^ dynes cm"^ from the VAL solar 
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1 /2 

models; 1 might be taken as the fractional area covered by 
evacuated magnetic fields. Strictly the coefficient l/yS, repre- 
sents the coronal base gas or non-force-free magnetic pressure 
extrapolated back to r which can be no larger than the average 
surface magnetic pressure. Thus the fractional area covered by 
evacuated magnetic fields at must be at least l/jS^^, or 1/71 
for Ps - 5000 and 1/32 for j3s - 1000, which corresponds 
to a minimum base coronal field strength of |Bcorl = 24 G for 
Ps = 5000 and |Bcorl ^ 54 G for ^ 1000, within the range 
of solar observations jLin et alJl2OO0l) . A lower transition re- 
gion occurs where the photospheric flux-tube areal coverage 
and coronal magnetic field strength are larger near active re- 
gions, consistent with the known tende ncy. 

The solar atmospheric model of iFontenla et alJ ( 119901) . 
which includes particle diffusion and conduction effects, lacks 
the Lya plateau and exhibits a much more abrupt transition 
region than what we obtain. While radiative, diffusive, and 
conductive losses must be largely balanced by MHS restor- 
ing flows in the corona, as we discuss in Sect. |2 loss mecha- 
nisms can influence the detailed shape of transition-region pro- 
files. Real radiative models coupled to more general magneto- 
hydro-steady equilibrium solutions need to be developed, but 
such work may be complicated by the intrinsic limitations 
and uncertainties in atmospheric modeling, especially asso- 
ciated with inhornogen eous magnetic structure tAvres 1981,: 
ICarlsson & Steinlll995l) . It is possible that general methods 
based upon transition-region emissivity profiles, which have 
been used to d etermine the nonradiative atmospheric heat- 
ing contribution dCraig & Browr]ll976l:lAnderson et aljfl 996). 
might be able to distinguish coronal heating by MHS restoring 
flows from other heating mechanisms. 

Figure|5]illustrates the solar atmospheric temperature func- 
tion T{r) for various values of ( with the coronal tem- 
peratures inferred from the white-light -intensity radial gra- 
dient taken from eclipse photographs jNewkirk et alJ fl97(il) 
and inferred from FeXIV 5303A line-width measurements 
j-Tarrett & von Kliiber 1958> . The limited resolution in r in the 
figure hides the transition region near r^. The coronal model 
for the spherical harmonic £ - \ seems to give the best over- 
all agreement with the measured radial profiles, but ( ^ 2.3 
better matches the base coronal temperature. However all of 
the model curves show a more rapid falloff than the mea- 
sured temperature profiles. Models containing an outward wind 
should give systematically higher temperatures with height in 
the corona. 



4. MHS restoring flows 

Under usual conditions, e.g. quasi-steady, the large-scale and 
long-term system average is represented by the equilibrium 
state. In the quasi-steady approximation, the system evolu- 
tion ca n be described by a sequence of nearby equilibria. Low 
ill977ll shows that the quasi-steady approximation is applicable 
when magnetic adjustments are fast compared to evolutionary 
processes and the magnetic field is dominant over the gas, as in 
the lower solar corona. 

A qualitative feature of MHD processes in a slowly evolv- 
ing dominant magnetic field where reconnections are fast is 




Fig. 5. Static solar corona T{r) for ^ =.5, 1, 1.5, 2, and 2.5 (from 
top to bottom) shown with temperature from white-light inten- 
sity in the quiet equator (filled circles) , in SW (cro sses), and 
in SE (open circles) streamers ( Newkirk e^ ^iJl^^T ^- and from 
line-w idth measurements (dashed line) (ilarrett & von Kliibeil 
Il958l) . 



an atmospheric restoring flow. Hydro-steady equilibria are al- 
ways self-restoring. In the absence of sufficient external gas or 
flow pressure, a time-dependent outward gas acceleration arises 
outside and along rigidly held non-force-free flux tubes in the 
overall MHD balance, as described in Appendix lEI Flows in 
a fixed magnetic field act generally to reestablish the equilib- 
rium atmosphere, heating the gas compressively. As long as 
the timescale for magnetic evolution is relatively long, such 
MHD processes can have a large-scale effect since the pres- 
sure perturbation that drives the flow propagates horizontally 
away from flux tubes at the sound speed. 

The available energy for outflows must be limited and the 
atmosphere has to eventually cool, distorting and geometrically 
stressing the coupled magnetic field, which contributes to mag- 
netic instability. The magnetic field can be perturbed in other 
ways too, for example by new flux coming up from the convec- 
tion zone below or by photospheric twisting and flux-tube mo- 
tions. Magnetic evolution is complicated since unstable equi- 
libria are generally possible where perturbations reduce the 
overall energy of the system. Also local deviations may per- 
sist within field allowances before equilibria can be resolved. 
However if the new magnetic equilibrium after fast reconnec- 
tion is mostly force-free containing rigidly held non-force-free 
flux tubes, outflows again must reestablish the equilibrium gas- 
pressure profile of a depleted cool atmosphere, applying the 
argument from AppendixlEl 

MHS restoring flows may provide an important coronal 
heating mechanism, which can exist consistently along with 
other heating processes. By offseting the total coronal en- 
ergy requirement, other heating sources should add stability to 
the large-scale magnetic field and thus longevity to the mag- 
netic evolution. The total coronal energy requirement is esti- 
mated conventionally using heat-flux scaling relations for loops 
and homogeneous re gions dRosner et alJI 19781 iHearn & KuinI 
198 1 ; Hammer 1982). If sufficient magnetic energy is not avail- 
able from evolutionary processes to offset losses, then pure 
force-free cool equilibria may be the only possibility. 

Localized heating in flux tubes, which is thought to be 
an important coronal heati ng process ( Vaiana & Rosner 19781 
lAirapetian & SmartJl995l) . may affect ambient MHS restoring 
flows. Heated flux tubes must exhibit field-line perturbations 
that stress the large-scale magnetic field and tend to produce 
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outflows in the surroundings at least on the average, which heat 
the ambient atmosphere above the normal equilibrium tem- 
perature. The ambient hydrostatic imbalance can have a rela- 
tively large scale of influence since it propagates horizontally 
away from its source much faster than thermal diffusion effects. 
Such post-reconnection hydrodynamic heating acting through- 
out the coronal volume may help overcome the often-cited dif- 
ficulty that reconnection events are so spatially localized and 
short-lived that only a larger event rate than what is directly 
inferred from obs ervations could accoun t for coronal heating 
(Hevvaerts & Prie slll984trHudsonll99ll) . 

Thus MHD processes may lead to refilling outflows, a 
'magnetic suction', which might be in evidence as the large- 
scale flows seen in the vicinity of coronal voids dWagner et alJ 
or be responsible for the supergranular-scale chro- 
mospheric upward velocities interpre ted from post-flare solar 
spec tral-line blue shifts ( Schmieder et al.lfl987t ICauzzi et all 
Il996 ). or appear as spicule eruptions in the chromosphere. The 
upward mass flux produced by solar spicule events is about 100 
times the total solar-wind mass flux, but represents only a small 
fraction of its total energy flux ("Pneu man & Ko pp 1978). 

The remarkable feature of an MHS corona is that the tem- 
perature is essentially a geometric parameter defined by the 
stellar surface gravity g,, radius r,, mean particle mass ill, and 
magnetic spherical harmonic f 



(23) 



2e + 4 

The observed average solar base coronal temperature of 1 .6 x 
10'' K corresponds to a spherical harmonic of { = 2.33 con- 
sistent with the large-scale solar magnetic field, which exhibits 
substantial power in the toroidal { - 2 and spherical harmonics 
^ = 3 or 4 as evidenced by the presence of active longitudes. 
With a general mixture of spherical harmonics, discrepant hori- 
zontal gas-pressure boundary conditions at flux tubes can arise, 
and a nonuniform large-scale surface temperature distribution 
is possible as a deviation from the pure MHS equilibrium state. 

Formulae for the average base coronal temperature like Eq. 
i23\ have been suggested akeady in the hypothesis of a 'geo- 
metric boundary condition' , relevant for different kinds of heat- 
ing models where the base pressure scale height is proportional 
to the stellar radius (Menzel 1968; Scudder 1992). The relation 
appears to gi ve reasonable coron al temperatures for a range of 
steUar types dWiUiams & MuUan>ll996) . 
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Appendix A: Local MHS equilibria 

There are well-known solutions to the classical static prob- 
lem describing a gas existing in conjunction with elongated 
magn etic fields without a perturbing gravitational field ( Parker 
Il979l Chapter 6). It is helpful to revisit the classical problem 
here in a general way. The classical MHS equilibrium equa- 
tions without gravity are the static equilibrium equation 



and Gauss' Law, Eq. @. 

Plane-parallel solutions can be developed based upon the 
general Cartesian vector function 



B 



da(x) da(x) db(x) db{x) 



dy 



dx 



dz ' dy 



(A.2) 



The X and z vector elements of B, da/dy and db/dy, are taken 
to be arbitrary functions, and the y element is written so that 
Gauss' Law is always satisfied. It is convenient to use derivative 
forms in the vector elements to avoid complicating integrals in 
the expression. 

Any solutions to Eq. iA.l\ must satisfy the two conditions 
B-Vp = and (VxByVp - 0, which are written with B from 
Eq. (|A3 



C„,a - Cy-b 
d 

dx 





dy 



'dz 



(A.3) 
(A.4) 



using the Unear commutation operators C.tv = If ^ ~ ^ ^t'^-' 
which contain p{x) as an implicit function. 

The two Eqs. (IA.3> and \AA\ can be used to write two of 
the functions in terms of the third, for example a{x) and b{x) 
in terms of p{x). The third relation, Eq. ( lA.U in the vector di- 
rection perpendicular to both B and VxB, gives a differential 
relation for the third quantity. The separate vanishing of ev- 
ery commutator term is a reduction that is satisfied with cer- 
tain general functional dependencies between the variables and 
leads to a well-known self-consistent form for the third rela- 
tion, the Grad-Shafranov Equation. Other possibilities are not 
considered here. 

Taking every commutator term to vanish in Eqs. ( IA.3> and 
(IA.4> . but without other reduction, gives a{x), b{x), and spatial 
derivatives of a{x) and b{x) as functions of p{x) alone, which 
is only possible with specific spatial forms for real functions, 
so other reductions are required in general. Taking d/dz = 
gives a consistent solution defined by the functional dependen- 
cies a = a(p) in Eq. (lAJt and f = f (p) in Eq. (|a3J. A 
similar solution is obtained by taking d/dx - 0. These two so- 
lutions or linear combinations are equivalent with rotation of 
the X - z axes around y. Taking z as the degenerate direction 
gives the usual 2D solution, ordinarily written with a as the 
implicit function 



B^i^^,-^^,BM^,y)) 



dy 



dx 



(A.5) 



Vp = — (VxB)xB, 

47r 



(A.1) 



The degenerate z is a direction of field elongation. The mag- 
netic potential function a(x,y) defines a 2D planform for the 
solutions, which determines the spatial variations in the lon- 
gitudinal field element B~ - B,{a(x,y)) and gas pressure 
p = p(a(x,y)) as arbitrary ID mappings. The magnetic field 
is everywhere perpendicular to the gradient of a(x,y), since 
B-Va(x,y) = 0. 

It is widely believed that this 2D solution is the only plane- 
parallel one. Using perturbative expansions, Parker ( 1979, Sect. 
14.2) shows that no nearby 3D solutions exist for bounded 
quantities in an infinite spatial domain. However 3D varia- 
tions do arise as large-scale deviations from plane parallel 
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jArendt & SchindleJli988l) . In the common shorthand, '2.5D' 
solutions refer to perturbed solutions in the problem with grav- 
ity, which admit large-scale variations in the degenerate direc- 
tion of the magnetic field. 

Appendix B: Properties of local MHS equilibria 

Substituting B from Eq. iA.5\ into Eq. JA.U gives the govern- 
ing equation for the implicit magnetic potential function a{x, y) 
for the static equilibrium problem without gravity 



V^a + — \4np{a) + 
\ da 



Vfl = 0. 



(B.l) 



Where Va 0, the potential function must satisfy the Grad- 
Shafranov (GS) Equation, written 



V^a(x,y) = -P'(a(x,y)), 



(B.2) 



defining the total pressure P(a) = A-Tip{d) + B,{a)^/2, where 
P'(a) = dP(a)/da . 

The total pressure P{a) and its component functions p{a) 
and B,{a) are all one dimensional and nonlinear, but must all be 
well-defined everywhere in the solution domain consistent with 
their forms at the boundaries. Strictly, disagreeing z boundary 
conditions are inconsistent with the z independence of the so- 
lutions; differences might produce small deviations in the solu- 
tions like twist or divergence or be a source of dynamical in- 
stability. It is popular to restrict consideration to entirely force- 
free solutions with p'(a) - 0; the restriction does not change 
the nature of the basic GS Eq. (IB.2> for a{x, y) but requires the 
specific current density / - j^B'Sa)B for an arbitrary B,{a), as 
is evident by expanding / = j^VxB using B from Eq. ( IA.5> 
with J oc B. 

The magnetic potential function a can be seen to be con- 
stant everywhere in the local plane of B and /. From Eq. 
dAJl . J-Vp{a) = 0. Then J-Vp{a) = {J-Va) p'(a) = 0. Thus 
/•Vfl = at least when p'{a) + 0, and when p'{ci) = 0, / || B 
so /'Vfl = anyway. We think of a as constant on ribbon-like 
current sheets, the solution surfaces being more curved along 
/ as it is defined by the derivatives of B. 

The GS Eq. ( IB.2t has the unusual feature that it contains 
the ID filter function P'(fl). The Fourier transform of the GS 
Equation in the linear case, with P'{d) - k^'^ft for a constant 
wavenumber, requires a transform function fl(fct, ky) that is zero 
everywhere in its 2D wavenumber domain k - (k^, ky) except 
on a thin annulus at |^| = ko, where arbitrary complex val- 
ues Hermitian in +k are allowed. Taking the Fourier transform 
back gives the potential function a(x, y) as a common Bessel- 
function radial kernel times azimuthal factor convolved with a 
spatial distribution of delta-function source points. Boundary 
conditions on a(x, y) constrain the azimuthal factor, leading to 
possible planar solutions defined by all power at one azimuth, 
or axisymmetric solutions with power uniformly distributed in 
azimuth, consistent with our visualization of current sheets and 
cylindrical flux tubes. Physical arguments show that nonlin- 
ear kernels are similarly c onstrained by boundary conditions 
llVaiiishtein & ParkeJl986l) . 

Random spatial distributions of a common GS kernel are 
nonlinear GS solutions too. We take the potential function 



a(x, y) to be the convolution of a distribution of sources D(x, y) 
with a common 2D kernel function A(x, y), a(x, y) = D(x, y) * 
A(x, y). The distribution is written D{x, y) - Yjj Cj6(x-Xj, y-yj) 
counting sources j of varying strength cj, where 6(x, y) denotes 
the 2D delta function. For a spatially incoherent or random 
distribution D of equal strength sources with all cj = 1, D 
raised to a power is D alone; for equal amplitude sources with 
Cj = +1, D raised to an odd power is D alone. For an ana- 
lytic nonlinear driver function P{a) the convolution factors out 
P'(fl) = P'(D*A) - D*P'(A), when the distribution inDis in- 
coherent, disjoint, and contains suitably restricted amplitudes. 
Then the GS Eq. JB.2I I reduces to the same GS Equation, but for 
the common kernel function A in place of a. Dark and bright 
threads appear to be isolated thin flux tubes aligned with the 
background magnetic field organized in current sheets sugges- 
tive of a bimodal amplitude distribution of a common flux-tube 
ker nel. 

|ParkeJ(ll972 ') considers MHS equilibria that contain a rela- 
tively strong constant background field. Expanding the Lorentz 
term in Eq. iA.l \ gives an alternate form for the MHS equation 



B' 



1 



(B.3) 



or with the superposition B(x) - B^ + B^(x) for a constant 
background field Be and a spatially varying field B^(x), we 
obtain 



(B.4) 



Taking the divergence, applying Gauss's Law V-B^ = 0, and 
using the property that the total pressure p + B^ji^n) is a 
bounded quantity in an infinite domain, we obtain 



{BeV)B-. = 0, 



(B.5) 



at least to first order in the small quantity Thus we 

have 



v(,.g).o. 



(B.6) 



with equality to first order in Eq. ( !B.5> proves to be 

accurate to all orders in the perturbative expansions, and the 
equilibrium Eqs. Q and are justified at least in the presence 
of a relatively strong constant background field, where the in- 
ternal flux-tube field is defined Bi = Be + B-.(- B). Additional 
discussion on Parker's perturbative expansion is contained in 
AppendixiCl 

The Parker theorem Eq. (IB.5> reflects the special features 
of the classical equilibria in Eq. ( IA.5> in the degenerate direc- 
tion of the magnetic field z- MHS equilibria allow only cer- 
tain alignments for the background field. A constant magnetic 
field can enter into the z vector element ^^(fl) without produc- 
ing other ramifications, whereas an added offset in the x or y 
elements of B requires adding a uniformly inclined plane to 
a{x,y). Such a plane limits the solutions and thus cannot be 
considered general, and even appears to be precluded by the 
boundedness of a(x,y): low- wavenumber components in the 
Fourier domain needed to represent the added plane are at odds 
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with the Fourier ring solutions in the linear case of the GS Eq. 
iB.2\ . and also contrary to the form of a superposition of spa- 
tially compact axisymmetric nonlinear kernels. 

Appendix C: Perturbative expansions 

The less-significant order « > equations in the perturbative 
series developed from Eqs. il3\ and il4\ are written for ^(r) = 
l;from the MHS Eq. O 

Vpi - i ((VxBo)xBi + (VxBi)xBo) = 

^(VrxBo)xBo, 
or for the general order n > 

V/3„-i;i:}=o(VxB,)xB„_,.= 

8^2:}:,;(VrxB,)xB„_,_i(l+2;^). 

From Eq. il4\ 



(C.l) 



2d 

or in general for « > 

V-B„ = -i-(Vr.B„_i) (l + 2(11 - 1)^ 
2d \ or 



(C.3) 



(C.4) 



The equations are linear in the variables p„ and B„ based upon 
order < n variables. The order-zero base functions po and Bq 
are 2D and constant in the degenerate magnetic field direction 
z; 3D eff'ects can be introduced with the cross-product terms or 
with general variation in the function h(r) with r not perpen- 
dicular to z- For consistency in Eqs. jl3> and J14> . Bo-Vpo - 0, 
which is the order-zero component of B-Vp - 0. 

A wide latitude is allowed for choosing perturbative quan- 
tities, and as long as they are relatively small almost every- 
where the standard procedure re mains va l id. Va n Ballegooijen 
( 1 19851 Appendix A) questions the |Parke3 d 1 972h choice of per- 
turbative quantity and equation separation procedure, but van 
Ballegooijen's assigning of approximate magnitudes to indi- 
vidual terms to separate his equations is not an algebraic pro- 
cedure, and his resulting equations have inherent contradictions 
as described bv >Parker (1987). Parker ( 197Z) does not ascribe a 
particular functional meaning to his kernel e, quoting from the 
beginning of his Sect. II: "Expand the field b and pressure p in 
ascending powers of some parameter e, which is of the order 
of e = where b is the small spatially varying magnetic 

field on the uniform background field B. Our choice for per- 
turbative kernel d/h(r) with d constant leads to a convenient 
separation of orders. 

The perturbative procedure is a functional separation 
method rooted in the natural independence of the basis func- 
tions e" irrespective of their real amplitudes. If the parameter e 
is small almost everywhere, the series solutions will be conver- 
gent. There could be locations where the amplitudes of compo- 
nents seem to be of diff'ering orders as van Ballegooijen con- 
tends. Where a normally small expansion parameter e actually 
becomes relatively large series expansions may be divergent, 
but such behavior characteristically defines the isolated singu- 
larities of differential equations as discussed in theorems on the 
Frobinius method for series solutions. 



Appendix D: Polytrope atmospliere 

A polytrope atmosphere is defined by the pressure relation 
p(r) oc p{rY for the polytrope adiabat F. Thus the ID poly- 
trope scale height h{r) - p(r) / (p{r)g{r)) can be written 



hs \p(fs)l g{r) gir)' 



(D.l) 



making reference to a base height rj taken to be the stellar sur- 
face, and introducing the scale-height function (p(r) for the nor- 
malized gas pressure. 

Substituting h(r) into the hydrostatic Eq. (|9j defines the 
scale-height function for a polytrope 



(C.2) <^(r) 



-i/r ^'f'i'') _ gjr) 



dr gshs 
. , , - F - 1 g{r) 

r = -, 

Or F gshs 



or 
d 



givmg 



F- 1 



Cgid 

J gshs 



dr 



and with g(r) - gsirs/r)^, we obtain 
F - 1 



(r) 



1 



F h. 



(-7) 



(D.2) 



(D.3) 



(D.4) 



(D.5) 



choosing the constant of integration so that <p(rs) - 1 . The part 
of the expression contained in brackets [ ] corresponds to the 
radial temperature variation, the thermal lapse rate. The poly- 
trope adiabat F ranges from 1 to 5/3 for an atmosphere ranging 
from lossy and isothermal to adiabatically stratified. For F = 1, 
the scale-height function <p(r) has a limiting pure exponential 
decay form. For F > 1, the thermal lapse rate and scale-height 
function have the unphysical feature that they become negative 
at some height, showing that a polytrope of constant F can exist 
only over a limited range in r. Such behavior is well known in 
stellar-interior polytrope models and leads to convection zones 
of finite extent. 

Appendix E: Hydrodynamics around fixed fields 

If the magnetic field is dominant over the gas, and either the 
field remains stationary, or the magnetic adjustment time is 
much shorter than the hydrodynamic adjustment time, then the 
Lorentz term in the MHD equation can be treated as constant 
in time and the hydrodynamics separate from the magnetody- 
namics. 

Fast magnetic adjustment is a remarkable feature of solar 
reconnection phenomena. S poradic cr ossings of coronal loops 
in low-energy flare events llLin et alJ suggest 'X'-type 

configurations, and the instability of the 'X' configuration is 
well known to form a new current sheet at near the Alfven 
velocity if a nonzero resistivity ca n be affected in the region 
('Duneev"l958'; 'Svrovatskii' 1981); such a resistivity is pre- 
dicted in some models (.Petschek & Thorne.l967;.PriesLl973) . 
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Taking a field strength of \B\ ^400 G as representative of chro- 
mospheric magnetic fields near neutral lines in solar active re- 
gions where flares commonly occur, we obtain an Alfven speed 
Ca = \B\/(4npy^^ ^ 10** cm s"' (with p = Nm = lO"'^ g 
using A' = lO'^ cm"^ and m = 10"^^ g). The Alfven speed stays 
fairly constant into the corona giving a collapse time of tb = 
hs/ca - 160 s for the propagation of the instability through a 
magnetic scale height lis - 2/zcoi- - rs/({ + 2) ^ 1.6 x 10'" cm, 
using a primary magnetic spherical harmonic of ( =2.33 in Eq. 
(EI}- 

The hydrodynamic response to akeady present or newly 
formed isolated thin non-force-free flux tubes can be obtained 
by perturbative analysis of the dynamical equations projected 
along the magnetic field. Taking the velocity to be aligned with 
the magnetic field and of lesser order than the mean quantities 
adds a single time-derivative term to the projected MHS Eq. 
Q, which is written with the continuity equation 



dui dp „ s ^ 

+ IT +P^(Vz- Vr) = 0, 

Ot OZ 



(E.l) 
(E.2) 



the equations are taken with the thermodynamic energy equa- 
tion p oc p^; z is the direction of the magnetic field; mi denotes 
the velocity amplitude; y is the ratio of specific heats, which 
ranges from 1 to 5/3 for lossy isothermal refilling flows to adi- 
abatic. We now choose to use the density p as a problem vari- 
able, rather than the temperature or scale height h. 

Perturbing the quantities backward in time, p - po-p\ and 
p - pQ-pi, from a final equilibrium state {po,po) leaves the hy- 
drostatic relation Eq. (|9} with po coupled to the magnetic field 
through the cross-field Eq. il 1\ . and three coupled perturbative 
equations for (p i , p i , m i ) 



dui dpi 
P^~^ " — -Hpi^cose, 
ot az 

dpi d 



(E.3) 
(E.4) 



and pi - yipol po)pi, introducing 6 as the angle between the 
magnetic field direction z and the radial r. With small positive 
pressure and density deficits pi and pi, a small positive ve- 
locity Ml is required in Eq. JE.3I I coiTesponding to an upward 
flow opposite the normal gravitational stratification. Refilling 
changes the gas-pressure profile, in effect heating the gas com- 
pressively as the dominant heating term in the thermodynamic 
energy equation. 

The unperturbed density po - po/(gh) varies on a scale 
similar to the gas pressure po, since the gravity g and tem- 
perature contained in h vary on much larger scales; thus the 
derivative is approximated dpo/dz po cos 0/h. The perturbed 
quantities must vary on the scale of the magnetic field 2h/ cos 0, 
giving d{p\,ui)ldz — > (pi,ui)cos6/(2h). Substituting for the 
time derivative d/dt — > 1/t„ lets us solve for the hydrodynamic 
timescale 



1 



4-h 



Tu = 



COS0 \ + 3y)g) 



1/2 



(E.5) 



which is proportional to the gravitational free-fall time through 
the scale height h. For a radial field = 0, mixed lossy condi- 
tions y = 1.2, and h - r J(2{+4), we obtain t„ =349 s for solar 
surface conditions with i - 2.33, which is just over twice the 
magnetic adjustment time Tg. The hydrodynamic adjustment 
time increases with more lossy conditions and smaller y, but 
overall depends only weakly on the thermodynamics. 
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